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Abstract — The rate optimization for wireless networks with 
low SNR is investigated. While the capacity in the limit of 
disappearing SNR is known to be linear for fading and non-fading 
channels, we study the problem of operating in low SNR wireless 
network with given node locations that use network coding over 
flows. The model we develop for low SNR Gaussian broadcast 
channel and multiple access channel respectively operates in a 
non-trivial feasible rate region. We show that the problem reduces 
to the optimization of total network power which can be casted 
as standard linear multi-commodity min-cost flow program with 
no inherent combinatorially difficult structure when network 
coding is used with non integer constraints (which is a reasonable 
assumption). This is essentially due to the linearity of the capacity 
with respect to vanishing SNR which helps avoid the effect of 
interference for the degraded broadcast channel and multiple 
access environment in consideration, respectively. We propose 
a fully decentralized Primal-Dual Subgradient Algorithm for 
achieving optimal rates on each subgraph (i.e. hyperarcs) of the 
network to support the set of traffic demands (multicast/unicast 
connections). 

Index Terms - Low SNR Gaussian broadcast channel, network 
coding, rate optimization, Primal-Dual Subgradient Method. 

I. INTRODUCTION 

Wideband fading channels have been studied since the early 
1960's. Kennedy showed that for the Rayleigh fading channel 
at the infinite bandwidth limit, the capacity is similar to the 
capacity of the infinite bandwidth AWGN channel with the 
same average received power [1, 2]. The robustness of this 
result in the case of with or without channel state information 
helps us model the low SNR wideband wireless networks in a 
general manner. It should be noted that when the band grows 
large, the transmitting power is shared among large degrees of 
freedom. This results in smaller SNR per degree of freedom. 
Using this as our underlying information-theoretic assumption 
to approximate the capacity over a link, we model the general 
traffic for this network and show that the linearity of capacity 
for disappearing SNR makes for the fundamental reason for 
simplicity in our model. Hence, we claim, it is possible to do 
networking over such a model with simplistic and essentially 
linear approach. 

In the context of wideband multipath fading relay channel, 
it was shown in [3] that in the non-coherent multipath fading 
relay channel, the same lower bound on the rate can be 
achieved as in the frequency division AWGN relay channel 
with the same received SNR, by using a peaky binning 



scheme. In this paper, we use a relaying scheme based on 
superposition coding. The rates achieved by peaky binning [3] 
are higher than the rates achieved by the relaying scheme that 
we propose. However, our relaying scheme has the advantage 
to extend easily to large networks due to hyperarc decom- 
posability. We would like to mention here that number of 
hyperarcs for Gaussian broadcast channel using superposition 
coding is equal to n for n receivers, instead of 2 n . 

The traffic model we use is quite general. It is divided into 
two classes: unicast and multicast (broadcast is considered as 
a special case of multicast), where each pair of source and 
receiver group in the network form a session for a particular 
class of traffic. The problem of successfully establishing mul- 
ticast connections in wireline or wireless networks has been 
long thought to be NP-Complete using arbitrary directed and 
undirected network models. With the advent of network coding 
(ref., [4], [5], [6]), and in turn breaking of the fluid model 
for data networks i.e. by performing coding over incoming 
packets, this approach has been able to intrinsically circumvent 
the combinatorial hardness of the multicast flow problem. 
It was also shown that establishing minimum cost multicast 
connections boils down to optimizing subgraph over coded 
packet networks [7]. 

In this paper, we consider a low SNR wireless network 
with Gaussian broadcast channel and MAC. The problem we 
consider is to optimize the rates for a given set of demands 
that needs to be met by this network. We show that this 
problem can be casted as a minimum cost multicommodity 
flow problem with intra-session network coding. 

This paper is organized as follows. Section II is composed 
of general problem formulation. In section III we propose a 
decentralized solution. We present our results in section IV 
and finally, we mention concluding remarks in section V. 

II. SET-UP AND PROBLEM FORMULATION 
In this section we introduce the hypergraph models for the 
low SNR Gaussian broadcast and multiple access channel. 

A. Low SNR physically degraded Gaussian broadcast chan- 
nel [14]. 

Consider a general wideband fading channel where the input 
waveform is x and the output waveform is y, the fading 
coefficient matrix is given by h and n is the additive white 




Fig. 1. (a): Two receiver physically degraded Gaussian broadcast channel 
with Zi ~ N(0, JVi) and z' 2 ~ AT(0, JV 2 - JVi). (b): Rate region for the 
channel in (a), dotted line denotes the flatness of the rate region curve in the 
limit of vanishing SNR with C'i and C'2 as max rates for each receiver respec- 
tively, (c): Decomposition into hyperarcs {(s, (<fi)), (s, (di, cfa))} with their 
common rates for the case in (a) with receivers di and di (corresponding to 
better and worse respectively). 

noise. The channel is given by: 

y = hx + n. (1) 

The capacity of the channel, for both Gaussian channels 
and fading channels increase sublinearly with the increase in 
signal to noise ratio (SNR) but in the low SNR regime the 
capacity in the limit is linear in SNR for fading and non-fading 
channels: 

C(SNR) = SNR + o(SNR)(nats/s/Hz). (2) 
Clearly at low SNR, the signal-to-noise ratio per degree of 
freedom (SNR) approaches zero in the limit [2, 8, 9]. For 
such a case the point to point capacity boils down to: 
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where h 2 = and D is the distance between the transmitter 
and receiver. Let us now look at the standard model of a single 
sender and two receivers with noise variances N\ and N 2 
respectively (ref. Fig 1(a)). The capacity region is given by: 

where C(x) = W(ln(l + x)), Ai + A 2 = 1, K > and P 
is the total power (ref. Fig. 1(b)), the transmission scheme is 
superposition coding [15]. 

The rate region defined in (4), when looked under the 
low SNR lens comes across as a rather simpler picture. For 
the power limited low SNR regime, the effect of the power 
allocated for the better receiver, as the contribution to the total 
noise experienced by the worse receiver is negligible (ref. Fig. 
1(b), for the rate region for low SNR in the limit). So, for the 
low SNR physically degraded Gaussian broadcast channel, the 
rate for the worst receiver can be approximated as 
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Generalizing the same idea for the case of a given source i with 
power Pi and n receiver nodes, where the receiver set J = 
(1, ...,n) can be broken into n subsets as J k = (1,2, ...,k) 
for k G [l,n] ordered in decreasing order of reliability. The 
rate region defined for each hyperarc (i,J k ) in the low SNR 
limit is given as 
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where, A^ < 1, which when combined appropriately gives 

fe=i 

the rate region of the set J k . The equation (7) comes from the 
fact that capacity is linear in the limit of disappearing SNR, 
where Lj for all i € [1, n] is the location of the node and a is 
the loss exponent. We formalize the above mentioned concepts 
and motivate our next definition. Let \ k Pi = Pij*, Vfc. 

Definition. 1: For a given sender i with total power Pi and 
a receiver set J = [l,K] in low SNR physically degraded 
Gaussian broadcast channel, the set J can be decomposed 
into K hyperarcs where each hyperarc is defined as the 
connection from the sender i to the receiver set J k = [l,k], 
where k C [1,-K - ]. The rate over each hyperarc is defined as 
r t j* = r U k = \\ Li _l'^ N2 , where, ^P iJk < Pi, Vfc e J 

k 

and the set J k ranges from best to worst receiver (ref. Fig 1 
(c)). 

B. Interference issues in multiple access at low SNR. 
Now, let's consider the case of multiple access where more 
than one node tries to access the channel at the given instance. 
Let there be U nodes in the system at an instance, and u C U 
of them are trying to access the channel at this instance, if 
node i e u intends to communicate with node j e U among 
others in u, the signal to interference and noise ratio (SINR, 
denoted as /xy) experienced at node j is given by: 
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Note that, since every node in u is interested only in a common 
receiver, we allocate the whole power of the node over this 
single hyperarc, so k = 1 and P ii7 i = Pi for every transmitter. 
But as we are operating in the low SNR regime, the intuition 
suggests that the effect of the interference should be negligible. 
We straightforwardly include it in our assumption, thus we 
define the rate (denoted with R) experienced at the receiver j 
as: , 
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The approximation (10) comes from the fact that the contri- 
bution of other signals being transmitted from other sources 
in the system with low SNR channel to the interference is 
negligible and the approximation (11) comes from the linearity 
of capacity in the limit of disappearing SNR (ref. Fig 2(a) and 
2(c)). In Fig. 2(b), we can see that the SNR curve approaches 




Fig. 2. (a): Two sender case for the low SNR multiple access channel, where 
Z ~ A^(0, N). (b): Rate region for case in (a), the dotted line denotes the 
respective SNR's fi\ and fi2 for two senders and the arrow shows that in 
limit of disappearing SNR, the SNR curve touches the capacity curve, (c): As 
the effect of interference is negligible, the case is (a) can be approximated as 
individual hyperarcs. 

the capacity curve in the limit, corroborating our assumption 
that the SNR equals capacity in the limit of disappearing SNR 
per degree of freedom. 

C. Low SNR network rate optimization. 
Let us represent the wireless network as a directed hypergraph 
Q = (A/", A), where Af is the set of nodes and A is the set of 
hyperarcs, where each hyperarc emanates from a node and 
a terminates at a group of nodes, which we also refer to as 
the broadcast group of the hyperarc. Note that we consider 
multicast in our multicommodity flow optimization model (as 
opposed to only unicast), thanks to network coding. 

It's important to note that the common rate associated with 
each hyperarc r itJ k, is the capacity of the hyperarc, because 
this is the rate that can be guarantied to all the receivers in this 
hyperarc. Also, r^k , is a nonnegative function of the transmit 
power Pijk of the hyperarc (i,J k ). Now that we consider a 
network with more than one sender, update of notations is 
required. For a sender i G Af, that is capable of reaching 
ki G (lj, .., Ki) nodes, where each ki G Af\i, the K hyperarcs 
are denoted by (i, J ki ), Vki G (lj, ..,-Kj). 

Imagine a set of traffic demands where m = 1 , . . . , M 
sessions need to be established, each with t m = 1, — , T m set 
of receivers, in a given wireless network that experiences low 
SNR and that is represented by the hypergraph Q = (Af,A). 
We know from the definition of hyperarc that a single node can 
lie on multiple hyperarcs, therefore, we need a way to carefully 
count the incoming information and outgoing information to 
apply the law of flow conservation to the hypergraph and 
finally be able to cast the problem as a flow optimization 
problem. For that, we define another graph Q' = (Af, A'), 
which is simply the equivalent directed graph of Q = (Af, A) 
with arcs instead of hyperarcs. This graph can be easily 
obtained by decomposing the hypergraph appropriately. Let's 
define the term (ref. [7] for detailed notation explanation): 

((ij k *)eA\J k i3i) 
which simply describes the way to add all the flow entering 



a node on all incoming hyperarcs, corresponding to the graph 
Q' = (Af, A'). Notice that x k { is not the same as r.,^ 
defined in the previous section, ar fJ * 4 can be interpreted as 
the flow between i and receiver I of the hyperarc J^, and 
it cannot exceed the common rate {r iJ k i ) associated to the 
hyperarc which is also the hyperarc capacity, for each I G ki. 
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ll-JVo - fiJ -^ iJki • Then, the mini- 
mum cost optimization problem for the low SNR network can 
be formulated as: 

minimize 
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where P» is given V i, x ™ k (m) G Frj k (m), and F:j k (m) a 
bounded polyhedron made of flow conservation constraints: 

Y z*^. (m)- Y x^.^rn) = Si(m), 

(j^\(i,j^)eA') ' (j^\(j^.i)eA>) 

Vi G 7V,Vt TO ,Vm 

(17) 

x t ™ k .(m)= ^ x u k i( m )> 

(j^£j k i\(iJ k *)eA) (18) 
V(i, Jj ') G A',Vt m ,Vm 
x>,(m) > 0,V(i, J^) G A',Vm,Vt m G [l,T m ]. (19) 

As opposed to standard multicommodity flow problem in 
which flows are simply added over a link, the constraint (13) 
in fact catches the essence of network coding by taking only 
the maximum among all the flows of a session (note that we 
only consider intra-session network coding). Since F*J k . (m) 
is the polyhedron formed by the laws of flow conservation, 
constraint (18) translates the flow conservation laws from the 
underlying directed graph A' to the hypergraph A (the wireless 
network) by adding the flows on all hyperarcs between node 
i and J-* i.e. flow in (i, J.*) G A' is the sum of all the flows 
on the hyperarcs (i, J fcl ), VJ fe * 3 J ki . 

As we can see, the above mentioned problem is a convex 
optimization problem. The only nonlinear constraint is (13), 
and could be readily replaced by the set of linear inequality 
constraints y lJ k i (m) > (x™ ki (m)), Vt m G [l,T m ]. The 
modified problem results in a standard linear multicommodity 
flow problem with linear objective and linear constraint set. 
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where a^™^ (to) e -^/jlj ( m )> ar *d P*Jki ( m ) is a bounded 
polyhedron made of flow conservation constraints. Note that 
we optimize the power over each hyperarc, to determine the 
optimal rates for each hyperarc that satisfies the network 
demands, we simply need to multiply the optimal power with 
7 iJ fc i . We will prefer to solve the problem by proposing a de- 
centralized algorithm for generally understood and appreciated 
reasons. 

III. DECENTRALIZED ALGORITHM 
For developing a decentralized solution for problem (B) we 
need to understand the structure of the primal problem first and 
transform it into a separable form. We know that the objective 
function is a linear and increasing in its domain and so are 
the constraints. 

Taking the Lagrangian dual of the problem (B) we get the 
dual optimization problem as: 

maximize^ £ QiJ k i + £ &Pi) (Q 

(i,J k i)eA idM 
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The dual problem is clearly hyperarc separable and could be 
solved in a decentralized manner. But the dual problem is not 
differentiable at all the points in the dual domain, this is due to 
the fact that there might not be a unique minimizer of q iJ k i for 
every dual point as the objective function is a minimum over 
sum of linear functions for fixed dual variables. To solve the 
dual problem (C), we need to solve its subproblem (25). The 
subproblem (25) (and the dual problem (C)) could be solved 
with a lot of techniques, [10, Chapters 8-10], [11-Chapters 5- 
6, 12-Chapters 6] using some subgradient based technique but 
they do not necessarily yield the primal solution (which is of 
our interest here). There are however, methods for recovering 
primal solutions from the dual optimizers. 



We will take a different technique than the above mentioned 
approaches but before lets look into some inter-dependence 
characteristics of the dual and primal problem structures. Sim- 
ply having convex primal problem in hand does not guarantee 
strong duality, but with some constraint qualifications we 
can assert that strong duality holds or not. One such simple 
constraint qualification technique is called Slater's condition. 

In our case it can be easily seen for constraints (13) (or 20) 
of problem (A) (or (B)), there exist a vector {x*J fc . (to)} for 
which the inequality can be strict. 

Let us represent the set of primal vectors as p = {x,y,z,P} £ 
Si where 5i is the feasible set for the primal problem, 
and similarly we can do it for the dual problem, d = 
{A, v, /i, £} £ 5*2- As we can see that the primal and dual 
optimal are equal (thanks to strong duality), we can express our 

problem in the standard saddle point form max min </>(p,d) = 

deS2 peSi 

min max0(p,d), where function (f> is the Lagrangian dual 

of the problem (B). This implies that for (C), we get the 
hyperarc separable saddle-point form 
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Now we are in the position where we can solve the problem, 
separable in hyperarcs using any saddle-point optimization 
method for non-smooth functions. For our problem set up, we 
propose a Primal-Dual Subgradient Algorithm by Nesterov for 
nonsmooth optimization [ref. 13]. Nesterov's method generates 
a subgradient scheme intelligently based on Dual-Averaging 
method which beats the lower case complexity bound for any 
black-box subgradient scheme. The algorithm works in both 
primal and dual spaces, generating a sequence of feasible 
points, and ultimately squeezing the duality gap to zero by 
finally approaching the optimal solution. A positive conse- 
quence of the Primal-Dual approach is that at each iteration we 
get a pair of points (p,d) which are primal and dual feasible, 
hence, we get the primal feasible solution with essentially no 
extra effort. As opposed to many subgradient type methods 
where there needs to be a method for primal recovery, specially 
for large and ill-posed problems. A. Primal-Dual Subgradient 

Algorithm. 

Since the dual function is hyperarc separable, we can optimize 
the power over each hyperarc separately and add each of the 
optimal solutions to construct the optimal solution of the dual 
problem (C), ultimately achieving the primal optimal solution 
for problem (B). The algorithm is as follows: 

1) Initialization: Set s = £ Q. Choose 9 > 0. 

2) Iteration (k > 0): 

. Compute g k = d(f>(p k , d k ). 

• Choose (Jfc > and set s k +i = s k + g k - 

• Choose 6 k+ i > 9 k 

( k 

Set y k+1 = ^-argm&x ^<J l (g(y l ),y l - y) 
xeQ 

where (g p ,gd) is the set of primal and dual subgradients and 
<r k ,s k and S k are aggregated sequence of points. 
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Fig. 3. Y-axis denotes 2 items, optimal primal costs computed using the 
infeasible path following method when applied directly to primal problem 

(B) and the optimal cost for the dual (primal optimal solution to (B)) when 
Primal-Dual Subgradient Algorithm is applied to the dual problem (C). 

IV. SIMULATIONS 

We now show the results of our simulations that support 
the claims of the algorithm presented. We solved the dual 
problem in a decentralized way by solving it for every hyperarc 
separately and then adding up the respective solutions to 
construct the dual optimal solution of the problem (C), when 
this is optimal, it is the primal optimal solution for problem 
(A) in our case. 

The setup consists of uniformly placed nodes on a chosen 
area of a x a m 2 , with given node locations. We start our 
simulations with smaller networks of only 4 nodes on a 10 x 
10 m 2 area with the area size increasing as the number of 
nodes in the network increase to keep the node density/area 
in a controlled range. Each node has a single hyperarc and 
it can communicate with all the nodes in the network, this 
is just a simple generalization of our case where a node can 
communicate with only a subset of total nodes in the network. 
For each network we randomly choose a set of m multicast 
sessions and T m set of receivers for each session respectively 
with the required rate demand associated with each session 
that need to be established, but making sure the the traffic 
demands are < the respective min-cut for each session to make 
the problem feasible. 

In Figure 3, we compare the optimal solution approxima- 
tions of the Primal-Dual Subgradient Method for problem 

(C) with the standard infeasible path following method for 
problem (£>). It can be seen that the our proposed algorithm 
gives close approximations of the primal solution of the 
problem (B). Note that the path following method is directly 
applied to the primal problem and the Primal-Dual subgradient 
method is applied to the dual problem, to compute the dual 
solution of the problem (C), which will be give us the close 
approximation to the primal solution of problem (A). 



and MAC. We do this by showing that rate optimization for 
the Low SNR physically degraded broadcast wireless network 
can be formulated as a standard linear multicommodity flow 
problem for optimizing power over each hyperarc using net- 
work coding. Our model is relieved from interference related 
issues, this is due to the fact that the capacity of the low 
SNR wideband channel is essentially linear in SNR per for 
vanishing SNR in the limit, which relieves the system from 
interference and related issues. Our model operates in the non- 
trivial feasible rate region that achieves capacity in the limit 
of disappearing SNR with appropriate encoding scheme. 

We use a primal-dual algorithm to construct a decentralized 
solution for solving the problem, which has apparent advan- 
tages for recovering the primal solution than standard projected 
subgradient methods. In the simulation results shown, we don't 
present the gains of routing using network coding over simple 
routing. But is already a vast literature establishing this fact. 

Finally, we believe that realizing low SNR networks is a 
worthwhile attempt as the linearity of capacity in the limit 
SNR — > provides a fundamental simplicity for networking 
to be done. Insights reveal interesting and promising work 
could be build up and blended with our simple model (e.g. 
mobility, reliability etc), which remains to be explored in this 
scenario. 
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V. CONCLUSION 
We develop an efficient optimization model that provides 
an achievable rate region for the low SNR broadcast channel 



